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Abstract. The form factor of a quantum graph is a function measuring cor- 
relations within the spectrum of the graph. It can be expressed as a double 
sum over the periodic orbits on the graph. We propose a scheme which allows 
fl ' one to evaluate the periodic orbit sum for a special family of graphs and thus 

to recover the expression for the form factor predicted by the Random Matrix 
Theory. The scheme, although producing the expected answer, undercounts 
f*^ ' orbits of a certain structure, raising doubts about an analogous summation 

recently proposed for quantum billiards. 



^ ' 1. Introduction 

in'. 

r\i , One of the central questions of quantum chaology is investigating which prop- 

er ' erties of the spectrum of a quantum system reflect the chaoticity of the underlying 

'nJ" , dynamics. It has been observed |H| that if the classical limit is chaotic, the local 

statistical properties of the spectrum of a generic quantum system resemble those 
of the spectrum of a large matrix from a suitable random matrix ensemble. The 
choice of the ensemble depends only on very general characteristics of the system, 
such as whether it is time-reversal invariant, and whether spin is present. Clari- 
fying the precise meaning of the term "generic" above and proving the conjecture 
of in] in any generality still remains one of the most exciting questions of quantum 
^ ' chaology. 

Quantum graphs make ideal models to study the questions of quantum chaol- 
%^ ' ogy, as they exhibit the same rich variety of behaviors and share the same set of 

analytical tools as the more complicated systems. However one can achieve deeper 
understanding on graphs while using methods which are simpler yet more rigorous. 
And due to direct analogies between graphs and more complicated systems like bil- 
liards and cavities, any advancement on graphs will provide insight to researchers in 
a much wider area. It is notable that some of the recent progress in understanding 
the role of periodic orbits in the spectral statistics was achieved on graphs first ^JH] 
and then extended to other systems J^. 

In the present paper we aim to exploit the analogy between the tools used in 
quantum chaology on billiards and graphs to analyze a recent breakthrough result 
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by Miiller et al [2I]- The authors of |23 investigate the form factor, one of the 
functions measuring the correlations within the spectrum of a classically chaotic 
quantum billiard. The advantage of the form factor lies in its convenient expansion 
in terms of pairs of periodic orbits of the classical system. This expansion is an 
application of the trace formulae (see |10| for a review) and has been a starting 
point for many investigations. The progress in understanding the role of periodic 
orbits in the universal behavior of the form factor of billiard systems is marked by 
such milestones as the diagonal approximation 6 , the first off-diagonal contribution 
|27II28| and, most recently, |2J| which announced a complete expansion of the form 
factor for the rescaled time t < 1. For quantum graphs, the corresponding steps 
have been taken in |18lll9im] (the diagonal approximation) and P) (the first off- 
diagonal contribution). This article aims to lay the groundwork needed for the 
complete expansion. It must be mentioned that the form-factor for graphs has also 
been successfully analyzed using non-perturbative supersymmetric methods |12ll3j . 

To introduce our results we need to briefly summarize the achievements of |21| . 
Using an expansion of the form factor as a double sum over the periodic orbits 
of the classical system the authors classify the periodic orbit pairs and perform 
the summation to recover, for r < 1, the result predicted by the Random Matrix 
Theory. The summation can be loosely divided into two parts: evaluating the 
contribution of pairs of orbits related through a fixed transformation (described 
as a "diagram") and summation of the resulting contributions over all possible 
transformations. While performing breathtaking (and mathematically rigorous) 
feats of combinatorics in the second part, the authors remain somewhat vague on 
the assumptions made about the ergodic nature of the periodic orbits in the first 
part. By analyzing the contributions of the diagrams on quantum graphs we aim 
to look for possible omissions in the derivation of 21 . 

We are pleased to report that, in the large, we are able to replicate and thus 
confirm the findings of 21 . In doing so, however, we are forced to omit a class of 
periodic orbits that do not fit well into the general summation. We believe that such 
orbits were also omitted (although less explicitly) in [2J|. These are the orbits with 
short stretches between the self-intersections. On one hand, if the conjecture of 
were to hold, such orbits should not contribute to the universal final answer: the 
short stretches are not universal. For exactly the same reason the task of dealing 
with such orbits is one of the hardest. A general method of proving that their 
contribution is negligible is badly needed for higher terms. 

The paper is organized as follows: in section [21 we discuss the definition of 
the quantum graphs, their spectra, the relevant statistics, the trace formula and 
the particular model considered. Then we review the past results and difficulties 
encountered in section |31 and proceed to describe our new summation method in 
section 01 Our results are compared with those of Miiller et al |23 in section O 
where we also speculate about the possibility of periodic orbit summation beyond 
the Heisenberg time t = 1. 

2. Spectral statistics of quantum graphs 

2.1. Quantum dynamics on graphs. Here we will describe a general ap- 
proach to quantization of a graph proposed by Schanz and Smilansky |24j . It shares 
many important features with another quantization procedure, used by several au- 
thors in the present volume (see, e.g., 53 )i which is based on finding a self-adjoint 
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extension of a second-order differential operator defined on the bonds of a graph, 
fn particular, in many important cases (such as the Neumann boundary conditions) 
the resulting quantum spectra coincide. 

We start with a graph G — {V,B) where V is a finite set of vertices (or nodes) 
and B is the set of directed bonds (or edges). Each bond has a start vertex and 
an end vertex, they are specified via functions s : B ^> V and e : yB — > V. The 
start and end vertices specify a bond uniquely: we do not allow multiple edges. We 
will, however, consider graphs with looping bonds, i.e. h such that s{b) = e{b). We 
demand that the set B be symmetric in the sense that 6 G ,8 iff there is a bond 
b E B such that s{b) = e(6) and e{b) = s(b). The bond b is called the reversal of 
b. A looping bond acts as its own reversal. The operation of reversal is reflexive: 

b = b. When we refer to the "non-directed bond b" we mean the couple of bonds, 
b and b. The number of vertices is denoted by 1/ = | V| and the number of directed 
bonds is i? = \B\. The vertices are usually marked by the integers starting from 0. 
The graphs we will be considering are metric, that is each bond b has a length 
denoted by Lb. Naturally, Lb = L-^. As a rule, we will be assuming that the different 
lengths are rationally independent, which means that the only choice of integers 
{kb}beB satisfying 

(1) Y^^bLb^O 

beB 

is the trivial one: fc;, = for all b E B. 

Now we define the quantum scattering map to be a i? x _B unitary matrix 

(2) U(fc)=D(fc)S. 

Here the matrix D is diagonal with the elements e*'^^'' , where Lb is length of the 
bond b. An element Sb2.bi of the unitary matrix S is zero unless the transition 
between 6i and 62 is possible according to the graph's geometry (i.e. e(6i) = 5(62))- 
The evolution is defined by ipn+i — U{k)ipn, where ip are _B-dimensional vectors 
of the wave amplitudes on the graph's bonds. The eigenvalues of the system are 
the wavenumbers k at which stationary solutions exist, i.e. U(/c) has an eigenvalue 
one. The quantum evolution is time-reversal (TR) invariant if the elements of the 
matrix S satisfy 

(3) 5*62,61 = Sb^M- 

The corresponding classical dynamics on the graph is defined as the Markov 
chain on the bonds of the graph with the transition matrix M^^-*, 

(4) Mb,M^\Ub,.bA^^\Sb,M\^- 

One of the most exciting questions of quantum chaology on graphs is to relate the 
statistical properties of the spectrum of the quantum graph to the ergodic properties 
of the corresponding Markov chain. 

2.2. Considered model, fn this paper we will only consider a particular 
family of graphs, complete Fourier graphs. In a complete graph, there is a bond 
between any two vertices, including a loop from each vertex to itself. In a graph 
with N vertices, B is equal to 7V^. The corresponding S matrix is defined by 

1 ■ 2x, 
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where j = s(bi) and m = 6(62) are integers between 1 and N. The distmguishing 
feature of a complete Fourier graph is the quick equihbration of the corresponding 
classical transition probability: the probability to go from a bond &i to a bond &2 
in t steps is equal to 1/B independently of 61 and 62 as soon as i > 2. 

2.3. Spectral statistics and the BGS conjecture. Consider an ordered 
sequence {fcij^i with mean density 

(6) -d^ hm Mn^^^J^ 

equal to 1. The most popular quantity for numerical study of the statistical prop- 
erties of the sequence {ki}^i is the nearest-neighbor spacing distribution, which is 
the weak limit (if it exists) 

1 ^ 

(7) pix) = lim ^ V ^(a; - (fc,+i - h)). 

i=l 

The two-point correlation function is a second order analogue oi p{x) and is defined 

by 

N N 

(8) i?2(x) = Jhn^-5]5]^(x-(fc,-fc,)). 

It is important to note that the normalization in front of the sum is N^^ and not 
N^'^, i.e. R2{x) is not a probability distribution. The most popular quantity for 
analytical investigations is the Fourier transform of R2{x), the form factor, 

/oo 
e2^"^(i?2(a;)-l)dx. 
-CXj 

It owes much of its popularity to its convenient representations via trace formulae. 

Example 1. Uncorrelated spectrum. Let {Ci}iZi be a sequence of inde- 
pendent exponentially distributed random variables with expectation 1. Define 
ki = X^fci 0' these are the jump times of a Poisson process [7j- O116 can show that 
R2{x) = 1-1- 5{x) for the sequence {ki}, i.e. there are no detectable correlations 
within this sequence. 

The Bohigas-Giannoni-Schmit (BGS) conjecture [H], if applied to quantum 
graphs, claims that the spectrum of quantum graphs is correlated in the same way 
as in the spectrum of large random matrices ^^EOI ■ In the case of TR invariant 
quantum graphs, the random matrix should be from Gaussian Orthogonal Ensem- 
ble (GOE) or Circular Orthogonal Ensemble^ (COE). In particular, the form factor 
of graphs, K^^\t) should converge to the form factor for the Gaussian Orthogonal 
Ensemble (GOE), 

(10) K{t) = 

as we increase the number of bonds B. Since often there is no unique way of 
increasing the size of the graphs, a separate conjecture describes the sequences that 




Loosely speaking, GOE corresponds to the spectrum of the Hamiltonian whereas COE 
corresponds to the spectrum of the S'-matrix. The spectral statistics are the same. 
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are expected to produce convergence to the RMT statistics. This conjecture, often 
called the "Tanner conjecture" was formulated in PJ (see also |12lli;^j ): a related 
criterion was proposed in [3f,1J. They can be summarized as requiring the mixing 
time of the Markov chain M'^^-' to grow sufficiently slowly with B. Our model, 
the complete Fourier graphs, clearly satisfies this requirement: the mixing time is 
independent of B. 

2.4. Trace formula. One of the main links connecting a quantum system to 
the underlying dynamics is a trace formula which expresses the eigenvalues of the 
quantum system through a sum over the periodic orbits of the classical system. 
For quantum graphs the trace formula was established by Roth in |22[l23j and 
rediscovered by Kottos and Smilansky in |18j . Before we can write it down, we 
need to give several definitions. Let Vn be the set of all sequences 

(11) p=[6i,62,...,6„], b,eB, n>2 

such that for any i = 1, . . . ,n the bond fe^+i follows the bond fe^: e(bi) — s{bi^i) 
(by bn+i we understand 61). Define the cyclic shift operator a on Vn by 

(12) a([6i,62,...,6„]) =[52,63,...,6„,6i]. 



Equivalence classes in Vn with respect to the shift a are called periodic orbits of 
period n. For a periodic orbit p = (5i, . . . , 6„) we define its length and amplitude 

by 

n n 

(13) Ip^^Lb- and ^p = JJ-Shi+i,!,. 

i=l i=l 

correspondingly. A periodic orbit p is a repetition if it can be represented as another 
orbit repeated r times. The maximum of all such r is called the repetition number 
and denoted rp. An orbit with 7'p = 1 is called prime. If n is a prime number, all 
orbits of period n are prime. 

The trace formula establishes a connection between the set of all eigenvalues 
{kn}^=i and the set of all periodic orbits V, 

(14) d{k) = V ^(fc - kn) ^^ + -y —Ap cos(fc/p), 

where L is the sum of all bond lengths of the graph. The right-hand side of H14|l 
is convergent in the sense of distributions. The trace formula is exact on quantum 
graphs and can be used to recover individual eigenvalues j5| as well as the whole 
density of states. 

2.5. Eigenvalue vs eigenphase statistics. The trace formula, (|14|l . pro- 
vides a convenient expression for the form factor of eigenvalues {fc„}^]^, 

where p (corresp. q) is a periodic orbit of length Lp (Lq) and with "stability" 
amplitude Ap {Aq). As signified by the last Kronecker delta, two orbits p and q 
will contribute to the double sum only if the lengths of the orbits coincide. 
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It is even more convenient to consider the form factor of eigenphases of the 
graph's scattering map. For a unitary operator, such as the matrix 1J(k) defined 
in ©J the form factor is defined at integer times n — 0,1, . . . by 

(16) X(f )(r) = l(|tr(U(fc))"p), 

where r is the scaled time r — n/B. In our case the averaging is performed 
either with respect to the lengths of individual bonds or with respect to k. Both 
lengths and k enter only via the matrix D and averaging produces equivalent results. 
Expanding the powers of the trace we obtain 

(17) tr(U(fc))" = ^-^pe 



p 



'^P 



Substituting this into the definition of the form factor and performing the averaging 
we arrive to 



KB)/ ^ 1 



(18) K'i^) = uLr^^p^q'^^p^^.- 

p,q 



B ^-^ rpTq 



It is an accepted fact^ that it is equivalent to study the spectral statistics of the 
sequence {ki} of eigenvalues and the statistics of the eigenphases of U(fc). Certainly, 
if we send all individual bond lengths to 1 in a fixed graph with real S, it is clear 
from expressions (|15|) and H18|) that Ks (t) will converge to Kjj ' (r) in the sense of 
distributions on any compact interval. 

We are interested in the limit of large graphs, B ^ oo, where Kfj '{t) is 
expected to converge to its limiting GOE form, 

^GOE(r) = 2T-Tlog(l + 2r) (0 < r < 1) 

(19) = 2t - 2t^ + 2t^ + 0(t'^) . 

When taking this limit we can ignore the orbits that are repetitions of other orbits, 
simplifying the formula for the form factor to 



(20) K^J'\r) = -J2ApA*^5,^,,^ 

p,q 

Equation (|20(l is the starting point for our derivation. 



3. Past results 

3.1. Diagonal approximation. To approach the task of evaluating the form 
factor one needs to classify possible pairs of orbits of equal length. The most simple 
such pair is p and q = p. Since we assumed that Lh = Lq;, we can also take q to 
be the reversal of p: q = (fo„, . . . , bi). It turns out |S] that such pairs provide the 
leading contributing term to the form factor, called the diagonal approximation. 
Evaluating it for graphs was studied in detail in [31] . Here we recap the derivation 
of the diagonal approximation for our model, the complete Fourier graph, for which 



although the author is unaware of a mathematical proof of such a statement, the reasons 
behind this correspondence are discussed, among other sources, in 1111191251 
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Figure 1. Notation around a self- intersection. The partner orbit 
is obtained by reversing the directions of the bonds number k and 
j and also reversing arc 2. 

the number of directed edges is TV^ and the modulus square of amplitude of any 
orbit of period n is N~". 



jy2 ^i^ V P p ' P p; ^2 

P P 



(21) 

n -T^ 2 „ n 

= ^7 > = 2^r = 2t. 

p 

Here we used the fact that there are asymptotically N'^/n orbits on a complete 
graph (we divide by n to take into account the shifts of the same orbit) and that 
the rescaled time r is n/B = n/N'^. 

3.2. The difficulties in evaluating further corrections. To appreciate 
the difficulties present when evaluating off-diagonal corrections, it is instructive to 
consider the simplest first correction in some detail. The first correction comes from 
the Sieber-Richter pairs |2ZIEHIj also called "figure of eight" orbits. We will say 
that an orbit p = (6i, . . . . 6„) has a self- intersection at vertex v if there are indices 
k and j such that k ^ j + 1 and s{bk) = e{bj) = v, see Fig. Q] Then the "partner" 
orbit 

(22) q= {bi, . . . ,bk^i,bj, . . . ,bk,bj+i, . . . ,bn) 

will have the same length as p (and it is easy to verify that the above q is a 
legitimate orbit). We do not require that k < j; the general rule is that we reverse 
all the bonds starting with bk and until we have reversed bj. 

To evaluate the contribution of such pairs we need to introduce more termi- 
nology, summarized in Fig.^ We will usually denote the intersection vertex by v. 
The sequence of bonds [6^+2, . . . , bk-2] we call arc 1; the sequence [bk+i, . . . , bj-i] 
is arc 2. The starting vertex of arc 1, 3(5^+2) will be denoted si, the end vertex 
e(6fe-2) — /i- Similarly for the second arc, s{bk+i) — S2 and e(6j_i) = /2. We 
define Ai to be the contribution to the amplitude made by traversing arc i, 

(23) Ai ^ Sbj^j^,bj+2-- ■ Sbk-2,bk~i^ A2 = Sbt,Mk+i- ■■ Sbj-i,bj- 

Using the above definition and the formula for Sb,b' we can write for Ap and ^q 



Ai X ——e^''^^ X A2X ——e^^'''^ 



(24) A„ = 



(25) A„ = Aix ——e^'f^^ x A-^ x —=i 
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(a) 



(b) 



Figure 2. Examples of orbits that would be miscounted in a naive 
summation. In orbit (a) the intersection vertex is not uniquely 
defined, choosing either v or S2 would result in the same partner 
orbit. In orbit (b) the partner orbit, obtained by reversal of the 
arc 2 coincides with the orbit itself, producing a pair which was 
counted in the diagonal approximation already. 



where A2 is defined by appropriately reversing the bonds in ^2- However, due to 
the reversal symmetry in the matrix S, Sb.t' = •S'p'bi we have A2 = A2. Therefore, 
the contribution of the pair (p, q) to the form factor sum simplifies to 



(26) 



^P^q 



i^in^2i 



1 



a(/lS2+/2Sl~/l/2-S2Sl)^ 



1 



Jf2-S2)(fl-Sl)i 



with the factorization in the exponent reminiscent of the expression for the action 
difference in billiard systems |30[I32| . 

It is now clear how one can sum over all "figure of eight" orbits. The sum will 
run over all possible intersection vertices, over possible choices for Si and fi and 
over possible configurations of arc 1 and 2 (including their length) . This approach, 
however, is prone to counting orbits that should not be counted. There are two 
basic corrections. 

The first correction comes from the orbits in which bk = bj, Fig.|2(a). In this 
case the orbit runs along itself and one can use S2 as another intersection point. 
It is clear, however, that this choice would produce the same partner orbit q. A 
solution proposed in |3] was to define the intersection point to be the leftmost vertex 
in such an "extended" intersection. In essence, this amounts to only counting the 
intersections with si 7^ /i. This condition was called a restriction. In the present 
paper we modify this approach. 

Another problem would arise if one of the arcs, say arc 2, was self-retracing. 
One of the simplest examples is when S2 = /2 and arc 2 contains just one bond, a 
loop from S2 to itself, Fig.|21(b). In this case, reversing arc 2 would give back the 
same orbit p and the pairing of p with itself was already counted in the diagonal 
approximation. Such orbits were termed exceptions in j2j; they were evaluated 
separately and subtracted from the previous result. We note that the previously 
imposed restriction si 7^ /i means that there can be no exceptions with arc 1 
being self-retracing. This shows that restrictions and exceptions cannot be treated 
independently. It should also be noted that both restrictions and exceptions must 
be treated carefully in order to get the correct result. 

To find higher order terms, orbits fitting more complicated intersection patterns 
have to be considered. Diagrams were used in U and to describe the relationship 
between the orbits p and q (i.e. which arcs get reversed). It proved to be hard to 
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Figure 3. The "figure of eight" orbits with intersection length 
0, 1 and 2. To ensure that the length is well defined, we impose 
restrictions on both sides of the intersection. The restrictions are 
drawn as double-ended arrows. 



find and reconcile all the exceptions^ and restrictions. We believe that the approach 
proposed in this paper helps to systematize the choice of restrictions in such a way 
that we avoid treating exceptions altogether. This approach was inspired by the 
recent work on billiard [21]. It has its own shortcoming which will be discussed. 
We believe that this shortcoming applies to |23 as well. 



4. A consistent approach to evaluating diagrams 

4.1. The second order diagram. We start by illustrating our approach on 
the simplest nontrivial diagram, the "figure of eight" . The trick is to introduce 
restrictions on both sides of the "extended" intersection and then sum over all 
possible intersection lengths, see Fig. |2| 

It turns out that we need to evaluate the first type of orbits on Fig.|31separately. 
The contribution of such orbits reads 

(27) i^2,o = ^E E ]^e-^/--)(/--)/-(l~J.,,,J(l-^.„,J, 

" Si,S2,fl,f2 

where we have implemented restrictions by using the Kronecker deltas. The term C 
counts the number of ways to complete the orbit by choosing the arc configurations. 
The two arcs are uniquely defined by the vertices they pass through. Between the 
two of them, they pass through n — 6 vertices (we have already chosen 6 vertices of 
the orbit: si, S2, /i, /2, and v twice). Thus there are P(n) = n — 5 ways to choose 
the lengths of the arcs and iV"^^ ways to choose the vertices themselves, leading 
to C = (n — 5)iV"^®. Going back to H27(l . we open one of the brackets to obtain 

ifon - — V V —e^'''(f^''''^^(f^'''^y^ (1-6 f) 

■" Si,S2,fl,f2 

(28) -^E E ^e-'/^-)(/--)/-^.,,,,(l-5.„,J. 

V Sl-S2jl-f2 

In both terms we perform the sum over /i. Since 'J2f e^^^^^^~^^'>^^/^ — Ss^j^ and 
^S2.f2 (1 ~ ^S2j2) = Oj the first sum is identically 0. In the second we substitute 



We extend the notion of exceptions to denote all orbits that could fit a particular diagram 
but that were already counted in another diagram's contribution 
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fi = Si throughout to get 

n^ v-^ v-^ C , ^ , n'^ ■,, , C 



(29) 



n^{n-5){N - 1) 



ATS 

The obtained result is out of order: since r = n/N'^ , we got an answer of the order 
t^N'^. Thus it is vital to evaluate the corrections from the extended intersections. 
It is easy to see that for orbits with an extended intersection, ApA* is indepen- 
dent of the orbit configuration and is equal to 7V~". For an intersection of length 
li, there are iV"~^~^'i(n — 5 — 2/i) possible arc configurations. There are also 7V'i+^ 
ways to choose the vertices denoted by the shaded circles on Fig. Inland N^{N— 1)^ 
ways to choose the vertices denoted by the empty circles. The total contribution is 

^ ]\[5+li ■ 

The contribution of all "figure of eight" orbits is thus 

n2(n-5)(iV-l) ^ n2(n-5-2Zi)(iV-l)2 



(30) 



TVS /l^ TV^+'i 

h=i 



-2— r- 

iV4 



Taking into account that r = n/N"^, we recover the correct second-order correction 
K2{t) = — 2t^. If we restrict the sum over Zi to the values up to (n — 5)/2, we get 
a correction term of the form 2n^/A^("+3^/^ which decays faster than exponentially 
as we send n and N to infinity. 

4.2. Evaluation of some third order diagrams. To gain some intuition 
for the general case we evaluate the contributions of two of the five diagrams con- 
tributing to the third order. The reader is referred to 4, for a detailed list of the 
third order diagrams and their "multiplicities" . The treatment of these diagrams 
was also performed for disordered systems in 29 , and for chaotic billiards in |15| . 

The first diagram we evaluate is the "sausage" orbits shown in Fig. ^ Here we 
need to require that all arcs (depicted by dashed lines) are at least one bond long. 
As noted by H. Schanz 26 this restriction excludes some potentially contributing 
orbits. Why they do not contribute to the form factor is a non-trivial question. 
In particular such orbits are not "rare" enough to guarantee a priori that their 
contribution is vanishingly small. We conducted a separate study to evaluate the 
missing orbits for this particular diagram and discovered that cancellations akin 
to those in H30() ensure they do not contribute. However, we are not aware of any 
general argument applicable to higher diagrams. 

In the intersections of length (first and second line in Fig.QJ we open up one 
of the restrictions, in the same way as in (|28|> . Another component necessary is the 



by that we mean an estimate based on the number of orbits and the modulus of an orbit 
pair's contribution, without taking into account the unitary cancellations 
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+ 2Z 



f^>^ 







+ 




Figure 4. The "sausage" orbits with various intersection lengths. 
The factor of 2 in the second hne accounts for the symmetric case 
h > 0, /i = 0. 




+ 




Figure 5 . One of the "clover leaf" diagrams contributing to the 
third order. In a triple intersection, such as above, a restriction 
prohibits coincidence of all 3 vertices at the same time. 



number of ways to complete the orbit by choosing the configurations of the 4 arcs. 
This is equal to P(n, /i, Za)^""^^"^''"^'", where 



(31) 



P{nMM) 



n-9-2li -2k 



is the number of possible choices for the lengths of the arcs — essentially the 
number of partitions of n into an ordered sum of 4 natural numbers. Now the total 
contribution of the diagram is 



(32) if3, = P(n,0,0) 



N^' 



7V2 



2 oc 



-2t72E^("''i'0) 



ii=i 



{N_rf 

]\[e+h 



^w. E "(".'..yir^ 



/l,i2 = l 



^3 n2(487V-32) „ , ^, ,,, 
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Next we evaluate a diagram with a triple intersection, shown in Fig. El The 
restrictions now mean that the three vertices cannot be all equal. In terms of 
Kronecker deltas it is implemented as (1 — Sy-^.y^Sy-^.y^). We again require that all 
arcs have non-zero bond lengths; this requirement excludes some orbits that ought 
to be counted. We must also note that in this case the restrictions do not fully 
rule out counting of orbits that ought to be excepted, for example the orbits with 
self-retracing arcs. The contributions of such orbits cannot be ruled out a priori — 
they occur often enough to contribute to the third order in r. However they do not 
contribute (as signified by the results of 01121) due to cancellations which should be 
studied separately. 

Turning a blind eye to these special orbits, we can evaluate the contributions 
of Fig. There are now P(ri, Zi)7V"~^~^'i arc configurations, where 

(33) PKM=(""V''^ 

Dealing with li — case separately we arrive to 

(34) '1=1 

- '^m+ N^N^-i) " -^^ [^ + <^(n ))■ 

4.3. The general formula. The big advantage of the new way of evaluating 
diagrams is that the restrictions we impose essentially decouple different intersec- 
tions from each other. This, in turn, allows us to write a general formula for a 
contribution of a diagram using only the information about the number of inter- 
sections and not how they are interconnected. 

Consider the "sausage" diagram again. To write its contribution in a more 
uniform fashion we introduce a special factor, defined for r > a > 0, 

1 if a = and r = 0, 

(35) M^ = <( if a = and r 7^ 0, 

otherwise. 

We can now rewrite the sums for K^a in the following form: 



Ka-l) 



(36) K,,^—yPin,r)Nf^ 



r=0 



r=l r=1 

where 

/n - 9 - 2r 
P(n,r)=(^ 3 

In the first summand of H36|l . r is a dummy variable since Mq — unless r = 0. In 
the second summand r corresponds to h in the second term of H32|l . Finally, in the 
third summand r — I1 + I2 and Mj counts the number of ways r can be represented 
in such a way. Now we can generalize further by denoting by V2 the number of 
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2-intersections in the diagram (in our case V2 — 2) and by introducing the notation 

v4 = 4w2 + 1 + 2r: 

a=0 ^ ^ r=a ^ ^ ^ 

Using the same transformation, the contribution K^h can be written as 
(38) ^3. = ^}:^(-l) (^J2.(^3 iJA^a(r) ^,,^^, , 

where ^3 = 1 and A = 6^3 + 3r + 1. From here we generahze in the foUowing 
fashion. Let, for a given diagram, the vector v — (v2, V3, . . .) describe the number 
of intersections of multiphcities 2, 3 etc. Let k be the maximal multiphcity in the 
diagram: Vj — for j > k. Then the contribution of this diagram is given by 



2 'U2 Vk k 

(39) ^(v)=^E--En(-i)"-- 

02=0 




J]M„,(rO(iV^-^-l) 



Vi-\-a,i 

•via,\ri) 1,JV - - ij 

i=2 



2=2 i=2 'i=2 






We studied sum l|5^ for several values of (w2,f3, . . .) and found that it gives 
the contribution derived in 21 , 



fc ^L-V+l 



(40) Z?(v) = (~l)^n^" (j.„y_i), (1 + 0(^^"')) 



i=2 

where 



E-< 

4=2 i=2 



Q = Jl r- and F 



At this point we can evoke the combinatorial results of ,21^ for the number of 
diagrams with the characteristic v. The summation of the leading order term 
of (|in|l over all possible diagrams then reproduces the predicted Random Matrix 
result, (Unj. 

5. Conclusions, comparisons with [21] and the outlook 

In the present article we outline a construction which allows us to recover ex- 
pression 1)4011 for the contribution of a general diagram to the form factor of a 
complete Fourier quantum graph. The summation over all diagrams, described in 
PT] would then recover the predicted Random Matrix result. However the final 
result, (|39|l . was obtained at the cost of ignoring the orbits which can be generally 
characterized as having short stretches between self-intersections. We believe that 
the same omission was used in |23 ; described in particular in the paragraph imme- 
diately following Eq.(16) (the discussion in Appendix D.2 is also relevant). While 
not claiming that the same applies to billiard systems, we can testify that on graphs 
the contribution of the orbits with very short stretches is absolutely divergent in the 
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limit B —f oo. This underlines the need to analyze the unitary cancellations which 
ensure that such orbits do not contribute to the final answer. The missing anal- 
ysis also holds the key to the possibility of extending the discussed summation to 
systems with slower mixing. 

Now that we have more understanding of how the universal contributions to the 
form factor arise from the sum of periodic orbits in the region r < 1, the question 
of why the expansion techniques break down beyond the Heisenberg time r = 1 
comes to the forefront. It seems, that if we are to recover the Random Matrix 
prediction for r > 1, the periodic orbit contributions would have to be classified in 
an entirely different way. Quite possibly, the summation needs to be based on the 
degeneracy classes — the sets of orbits of the same length — rather than on the 
diagram partition which is in a way transversal to the degeneracy class partition. 
Such degeneracy-based summation has been performed before for a special type of 
graphs, the Neumann star graphs [2] (see also (TCi; for an overview of the results). 
The radius of convergence of the expansion was found to be finite but, if re-summed 
^, the result was valid for all values of r. Important input for the question of going 
beyond r = 1 can also be obtained by comparing the periodic orbit expansions with 
the non-perturbative methods of |12lll8j which are valid for all r. 
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